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ABSTRACT 


This document develops transfer matrices which permit analysis and computer evaluation of 
(xrtain series-type microwave antenna circuits associated with an L-Band microwave radiometer 
(LBMR) under investigation at Goddard Space Flight Center. This radiometer is one of several 
diverse instrument designs to be used fcr the determination of soil moisture, sea state, salinity, 
and temperature data. Four-port matrix notation is used throughout for the evaluation of LBMR 
circuits with mismatched couplers and lossy transmission lines. M^'.rix parameters in examples 
are predicated on an impedance analysis and an assumption of an array aperture distribution. 

The notation presented here is easily adapted to longer and more varied chains of matrices, and 
to matrices of larger dimension. 
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TRANSFER MATRICES FOR SERIES-TYPE MICROWAVE 
ANTENNA CIRCUITS 


INTRODUCTION 

This document records some of the preliminary work which has been done to analyze the 
series-type circuit associated with the L-Band Microwave Radiometer (LBMR) design at Goddard 
Space Flight Center. The matrix formulation presented here is not restricted to the LBMK, and 
may be utilized in a general way. Furthermore, additional parallel paths may be obtained by 
increasing the dimensions of the matrices. A relevant example may be found in Ref. 1. 

The transfer matrices in this document are defined over a field of complex numbers. They 
are a direct consequence of the definition of a scattering matrix of some microwave component or 
junction. See Ref. 2. The latter may be postulated or may be obtained by means of a direct 
physical measurement on a “black-box." 

A general or bilateral transfer matrix may be regarded, after appropriate rearrangement and 
partitioning, as being composed of two null matrices and two unilateral (“transmission" and 
“reception") matrices of equal dimension. This paititiouing has been exploited throughout the 
present document; only unilateral “transmission" matrices are used. In this manner 7S-percent 
of the matrix -element bookkeeping is eliminated. The notational conventions, conversion from 
scattering to transfer matrix, and dimensional reduction are detailed in the beginning of this docu- 
ment. The development is then carried to the design stage for the overall problem. 

TRANSFER MATRIX CONVENTIONS AND REDUCTION 

The conventions used in this document are those of Ref. 1, and are made clear by means of 
a simple example: two unequal lossless line-lengths in parallel. 
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An output column matrix (g) results from the product of a transfer matrix (T] with an input 
column matrix (B). 
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Figure 1. Transfer-Matrix Conventions. 
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it follows that the reduced matrix for transmission from (L) to (R) is 
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( 2 ) 


MATCHED DIRECTIONAL COUPLER 

The scattering matrix (S) of the matched directional coupler of Ref. 2, p. 301, is rewritten, 


before reduction, as 
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(3) 


in view ol the port designations of the previous discussion. By inspection, the reduced transfer 


matnx for “transmission" is 
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Cl 
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(4) 


and represents the foui-port junction shown below. 
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Figure 2. Matched Directional Coupler. 

Since the matched*coupler matrix is both symmetrical and unitary, 

C,2+C2^*1 . (5) 

See Ref. 2, p. 301, and Appendix A of this document. Other [T] forms may be written for the 
matched coupler as there is arbitrariness in the phases of the signals at the output ports, depending 
on the choice of output reference planes. 

LOSSY TRANSMISSION LINE 

The transfer matrix for the lossy transmission line is predicated on equations found in Ref. 3, 
p. 13. Signal flow is assumed to be from left (L) to right (R) in Fig. 3, below. 
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Figure 3. Lossy Transmission Line. 


where 


and 


V _ ^ A e'T'* + B e'y* 

- Z_ ■ , 

^ A - B e"^* 


z = -fi 


7 = a-*-j/3 . 

A voltage reflection coefficient at a generic point on the tine is now defined, using 
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iac 


= Ae'>® 


( 6 ) 

(7) 

( 8 ) 


and 




( 9 ) 


as 


• \ 

At the special point K - 0, this reduces to 


Vine. A 


^ B • ^0 

r_ = — * .1 

° A Z. + Z, 


L *'0 

The line voltages at fi # 0 and 8 = 0 may then be wr tten as 

V - A t \ c-^8 - - a cosh t 8 ^ sinh 7O 

\^L ^ ^ 0 / '*’^' 0 ^ 


and 


Vl = A cTfi + B e^® = A + B A + A 
respectively. From these, after some algebra, the useful form 

[cosli 76 + 




Vl = 


(z[)“"h^' V„ 


( 10 ) 


(r.) 


( 12 ) 


(13) 


(14) 


IS obtained. 

An alternative to equation ( 14) is obtained more easily by writing the line voltages at 8 # 0 
and 8 <= 0 as 

Vin = Ae>« + r„ Ae-T« (15) 

and 

= A + r^A»(l +r„)A , (I6> 


respectively. From these, the useful fonn 


V. = 


(nr„)e-T'»V„ 


*■ (i»r„t-=Te) 


(17) 


4 


is obtained. Equation ( 1 7) is now incorporated into a unUaterai transfer matrix for a pair of lossy 
lines since two lines will be associated with the coupler output and since the entire series-type 
microwave circuit for the LBMR is cast in four-port notation. 

K is noted that both (V^„) and (V^^) above ate, in general, standing-wave voltages on the 
transmission line. When Pj, “ 0, only a single (traveling) wave is present. 

0 

(I +r„, 

where subscripts ( 1 ) and (2) refer to the uppc; and lower lines, respectively, and the (Z^ ) used 
previously goes over to (Z^^ j ) and ) for ports (3) and (4). 





MISMATCHED DIRECTIONAL COUPLER 

A unilateral or "transmission” transfer inalriv ir the mismatched directional coupler is now 
predicated on the reflection coefficients at coupler output ports (3) and (4). The conditions of 
Figure 5 are assumed; Zj^ j # Zq, and Z^^ ^ Z^ in general. The unitary transfer matrix of the 
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matched coupler is retained in the development. It is also assumed that there will be only a single 
input h| ^ 0 at port ( 1 ). and a termination ^ at port (2) of the coupler. 



Figure S. Mismatched Directional Coupler. 


The voltages across (Z| and (Z^ are the result of an incident and reflected wave super* 
position identical to that at tlie load-end of a mismatched transmission line i'i = 0). lu]uation ( 16) 
suggests that an auxilliary matrix should follow the classical unitary transfer nratrix of the coupler, 
and that it should have the form 

(M-r,,) 0 

IT) = (20) 

to effect a conversion from traveling to standing-wave voltage. This notation provides for indepen- 
dent complex-valued relleclion coefficients at ports (3) and (4) which will in general differ frmn 
(F^ j ) and ( (^ 2 ^ associated with a pair of transmission lines to the right ( K) of a coupler. Uqualion 
(20) may also be obtained by a fundamental argument based on energy conset/ation and may be 
displayed graphically on a Smith impedance chart. 

In the event that Oj = Oj = 0 , 

ir„i*iFo,i (21) 

and 

ir,2i = iro2i . (22) 

but complex equality is obtained only when K, * fij equals some integer multiple of a half-wave 
length. 

Since matrix products are not commutative the order of matrices for the case of mismatch 
must be as shown. 
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( 23 ) 


That is 



0 




jc. 






{T,,.l=lT„HT„.j]...[T,l (24) 

when the matrix subscripts increase in going from (L) to (R) in traversing a circuit. 

It is noteu that (23) serves for the unilateral (L) to (R) “transmission” under hj = 0 whereas 
matrix (4) allowed hj # 0. A brief discussion of the consequences of multiplying the two square 
matrices of (23) with a view toward obtaining an equivalent matrix under hj = 0 and an inquiry as 
to wliether or not the result is unitary, may be found in Appendix A. 

it is also noted that (hj) above is not a standing-wave voltage, but (g^) and (gj) may be 
standing-wave voltages. 


CONVERSION MATRIX 


Equations (15) and (16) show that (Vj^^) and (V^) are standing waves. Equation (23) shows 
that (hj ) is a single incident wave. If a self-consistent chain matrix is to be developed, an auxilliary 
matrix will be required to convert (Vj^) across the input impedance (Zj^^) to a coupler to an (hj ) 
implicitly across an impedance (Z„). Since (Z.^^) for the coupler at port (1) gives rise to a reilection 
coefficient. 




Z - Z, 

in 0 

Z + 7 


( 25 ) 


the inverse of the process leading to equation (20) results in a matrix of the form 

, Si')*’ 0 

IT) = 


(1 + S, )* 
0 


0 


( 26 ) 


which effects a conversion from standing to traveling wave voltage. It remains to develop an eiiua- 
tion for (Z^^). 

In analogy with transmission line theory, three waves superimpose at port (1) so that 
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(27) 


= Zo 


A + B, 


A-Bj-Bj l-C^p +c^r ** 

• ^ ‘ '^1 ' rl ^ *-2 * t 2 

The power into obviously equals the total power appearing at points (1), (2), (3), and (4) 

upon taking all internal coupler reflections into account. On the left side of the coupler, the 

relati ve power may be assessed using the equations 


IC'i^ r„ -C,2 r,2|2 . port(l) (28A) 

ye, (\i* ir,| . porto i:«U) 

On the right side of the coupler re^l^tiv^ power may be assessed us>ng the equations 


III + C|* 1 (R^ Z, , port(3) (2‘)A) 

111+ r,rr f\’l tR, Zj ^/|Z, j’) . port (4) . C‘)B) 

linergy conservation is easily venfied by the preceding, and absolute power division may also be 
determined. 


atti:ni!ation matrix 

four-port coupler with one input port terminated in its characteristic impedance is 
sometunes regarded as a three-port junction, leading to certain conformability problems regarding 
matrix multiplication. One convenient means of removing one of the coupler outputs on a chain- 
matrix calculation is to constnict an attenuation matrix, 



( 30 ) 
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which may be regarded as a specialization of a general attenuation matrix 


with Pj j = 1 , and 2 “ f'l 2 “ 1 “ ® ■ 


T = 

r,, 

I'll 



^22 

= 0 . 


•4 


( 31 ) 



Four-port notation may then be used throughout the analysis. The attenuation matrix may be 
introduced immediately after the coupler or some other component, such as a length of trans- 
mission line, may precede it. In any event, the order of equation (24) must be preserved since the 
matrix products are, in general, non-commutative. 

PRACTICAL CIRCUITS (UMR) 

A restricted example is now given to illustrate the type of circuit associated with an LBMR 
antenna array. Only four radiators and three couplers are considered. See Fig. 8. All transmission 
lines are assumed to be lossy, and the loss parameter (a) may be different for each line of arbitrary 
length. Both antenna impedances and loss parameters may lie functions of frequency. In the event 
that antenna impedances depart from the characteristic impedance (Z^,) of the lines, the loss 
parameter is redefmed, in the presence of standing waves, as 

(1 + |Pol*) 


Of = Of 


( 32 ) 


(i-ir,,i‘) 

where = a in Ref. 3, p. 3 1. Also see Appendix B of this document for a relationship between 
decibels and nepers when In equation (32), the magnitude of the reflection coi-fficient is 

taken to be 


IPol = 


\ ~ ^0 

h +^0 


(33) 
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It is tacitly assumed here tliat all lines have the same characteristic impedance (Zg) although 
the restriction may be lifted it necessary. A common notation is employed for all couplers, lines, 
etc. The values associated with tne upper and lower feed system are designated ( I ) and (2), 
respectively. It follows tliat (a, ) and (Cj ) of one matrix may therefore Iw different from (a, ) and 
(t| i of some other matrix. The circuit is now represented hy (n) matrices and the priKluct 


r ||„| Ii„ ,1 n„ ,.I on 

Here n = 1 5. Output may be obtained at any circuit interface by truncating the chain of matrices at 
the appropriate matrix 



Figure 8. Representative LBMR Circuit. 


The individual matrices are now written explicitly 


(1 e''''''' 
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(35) 
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( 37 ) 
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ITJ = 


(i+r„) 
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0 




ITJ = 


(l+r„)e''>'' 




~(T,„1.|T„I 


•72 «2 


|TJ = 


1 
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(I+ro2)e 




0 


'-IT,,! . 


(38) 


(39) 


(40) 


It is noted that (Tj 1 is similar (~) to (T^ 1 , (T, j 1 . etc., and not identical since the matrix values are. 
in general, distinct. 


PARAMETERS OF CIRCUITS 

A procedure is now outlined for obtaining the parameters of a circuit such as the one shown in 
Fig. 8. The procedure begins at the extreme right (R) of the network, and works progressively 
toward the left (L), with an impedance analysis. A discussion of pammeter determination unaviod- 
ably infringes on system design, which is beyond the scope of this document. Certain simplifying 
assumptions are made, however, it will be shown that these are realistic, not particularly restrictive, 
and may be lifted if desired. Further, it appears that a multi-stage procedure is r'xpedient for 
sopliislicated problems. 

Initially, assume the following: 

( 1 ) All line characteristic impedances are equal to the value (Z^), which is taken to be real. 

See Appendix C of this document. 

(2) All load impedances are also equal to (Zq ). 


II 



(3) The loss parameters (a) are equal to zero. 

(4) The current distribution is known for the array radiators. 

(5) The phase distribution is known for the array radiators, and is taken to be a constant 
here for convenience. 

(6) Coupler values (C| ) and (C,) will be predicated on the current distribution muttnitude, 
ignoring line loss (set a = 0), standing wave loss (set - 0), and rellection loss (set 

r, = 0). 

Under the preceding assumptions, the individual matrices simplify to 


IT, 


IT,1 = 


IT,) = 


(TJ = 


ITJ = 
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0 
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1 
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0 • 

e ^ 
0 ' 
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-IT^l.lT,^) 


-ITJ.IT.j) 


|T,1,IT,J 


-IT.oMT,,) 


1T,,1 


(41) 


(42) 


(43) 


(44) 


(45) 


(46) 
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Suppose that a current distribution (K^) is desired for (i) radiators, and that (i-1) couplers 
are to be used. See Figure 9. 



Figure 9. Power Distribution of Radiators. 

Linder impedance-match conditions (Z^ , - j = . . . = . = Z^), the radiated power of each 

antenna is proportional to the square of the input current, 

P. a|K,l^ = K.2 (41) 

A rule can then be improvised to determine the coupler values for the initial set of assumptions. 


“E 

(Coupler 1) 

(42) 

2 

a K,‘ 

(Coupler 1) 

(43) 

1 

a V K 2 

/S ' 

(Coupler 2) 

(44) 

3 

a Kj" 

(Coupler 2) 

(45) 

a K,^ 

(Coupler i - 1) 

(46) 


(Coupler i - 1) 

(47) 


As an example, when the current distribution is uniform, and 

Kj* = Co lant = 1 , 

and i = 4 as in Fig. 8. the coupler parameters are obtained via 

C|^=3C^^ (Coupler I ) 


(4X) 


(49) 
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(Coupler 2) 


(50) 


and 


C,2 *2C2^ 


C 2 = c ^ 
'“1 '^2 


C,2 +C2^ * 1 


(Coupler 3) 


(51) 


(5) 


Since a = 0, matrices IT, 1 , (Tj 1 , [Tj q 1 , ITj j 1 behave as simple phase shifters and the line lengths 
(K) arc selected to obtain a zero phase gradient (or any other) across the array aperture by regarding 
the complex operator of the coupler matrix as a phase advance: 

i »/2 


j=e' 


(52) 


Next, assume the following; 

( 1 ) All line impedances are equal to taken to be real. 

(2) Load impedances 1 > 2 > - * - ^Ai complex values not necessarily related to (Zg). 

(3) The loss parameters (o) for each line are not equal to zero. For convenience assume (o) 
is the same for all lines. 

(4) The square root ol tlie radiated power will be regarded in the subsequent array analysis, 
and will replace the antenna current distribution used earlic'. See Ref. 4, p. 142, which 
represents the electrical far-field as 


H = 


1 » . „■* Si ' 

lad . C 




-I r 


giO,0) . 


(53) 


(5) The design power distribution will be retained, but will be modified to account for ohmic 
or joule power loss due to standing-waves on transmission lines with a # 0. 

(6) The phase distribution is known for the radiators and is taken to be a constant for 
convenience. 

( 7 ) Coupler values will be predicated on amended power distribution which anticipates line 
loss (Oj), but not on reflection loss (Pj) In practical situations where the Z^ , may vary 
with frequency, rcllection loss is ordinarily just accepted since it would be impossible to 
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design couplers to compensate over a band. In fact, the assumed coupler properties of 
equation (4) arc a narrow-band approximation at best. Widc-band and single-frequency 
match are beyond the scope of this document. 

Given the second set of assumptions, it is then possible to begin with the antenna impedances 
(Z^^) and (Z^j), evaluate the reflection coefficients (Pj,,) and (rQj)and mismatch line 
parameters (a^ , ) and equations (1 1), (33), and (32). From equation (19) the imped- 

ance values at ports (4) and (3) can be determined. These, in turn, lead to evaluation of reflection 
coefficients (I\, ) and (P^j). From equations (27) and (25) the reflection coefficient (1\', ) at 
port ( 1 ) is easily obtained using the (Zj„) resulting from trple-wave superposition. Since 

r.', (541 

for the transmission line matrix ITjq) at coupler port (>), and since parameters7, , K,, 
the line matrix are already known, all parameters of (Tj j I throu^t (T| q ) are known with the 
exception of those in (Tj j 1 . 

From the relation 

(P +6P.) 

dB = 10 1og — — = 8.686 a 8 (55) 

10 (s) 


the power loss due to standing waves may be assessed, and the design power distribution value tor 

the i*'* feeder amended. The coefficients of the rightmost coupler (3), may then be set using logic 

similar to that leading to equations (46) and (47) once the losses of the lines (1 ) and (2) of matrix 

(Tj j I have been determined. Since equation (19) may be used again to determine impedances 

(Z|^ j ) and (Zj ^ ), with their associated relfection coefficients (P,, ) and (P, 2 ^ 

a’spectively, of the next coupler (2), the indicated process may be continued until matrix |T| 1 is 

reached. 

The standing-wave voltage (V^i^) at the left end of the transmission line described by matrix 
IT, 1 is given l>y 


V 


in 


( 


^0 1 I / 


*NI I 


(5t.) 
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where {: f ) is the input impedance to the entire network as seen to the riglit of the generator 
terminals. See Fig. 10. 



Tire driving function for the entire problem is, therefore, a column matrix 


0 


which operates on matrix (T| ] in the non-commutative sequence 



(57) 


(58) 


Returning to the discussion of coupler coefficients, suppose that an amended power dis- 
tnbution such as the one depicted by Fig. 1 1 has been developed by means of equation (55) and 
the original design distribution. 



Figure 1 1 . Amended Power Distribution. 


16 



The power (P| fiP^) is associated with the i"* feeder, and the increment (6Pj) is intended to offset 
ohmic feeder loss in the presen(% of standing waves. Due to uncompensated reflection loss, the 
actual power arriving at the P** radiator may be less than or equal to the design distribution value 
(P|). Reflection losses may be at ports ( 1 ) and (2) of the couplers, depending on values ^ ) and 
(Zl4)- 

Equations analogous to (46) and (47) can then be written as 

2 

R, a P; + 6P. (59) 

Lj I I 

2 

Rl^ a P„ +6P,., (60) 

The system of equations is solvable, using equation (S), and the summing process is carried out as 
before. 

CONCLUSION 

This document presents an approach to retaining transfer or chain-matrix methods in the 
presence of standing waves of voltage. An impedance analysis preceded the multiplication process 
to establish the matrix parameten. Current was not carried explicitly in the development. 

The analysis of the mismatched coupler was predicated on the notion that its classical 
scattering-matrix properties are invariant with respect to impedance mismatch. On this premise, 
reflected signals were described in terms of complex reflection coefficients and superimposed in 
analogy with standard transmission line theory. Energy conservation was verified for the coupler. 
The analysis was restricted to a single input to the coupler. It was found that for the general mis- 
match case, two waves superimposed at each output port and three waves superimposed at the 
input port. 

A brief discussion addressed the problem of determining the coupler coefficients for an 
assumed antenna array power distribution. The radiated design power from each antenna was 


c, (1+r,.) 


c,(i+r,p 
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used in the calculations since antenna input impedance was not necessarily identical among the 
array elements. Line losses were compensated in the determination of the coupler coefTicients. 
Impedance mismatch losses over a frequency band wen. not compensated. 
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APPENDIX A 


UNITARY MATRICES 

The term “unitary matrix” appears from time to time in the discussion on scattering and trans- 
fer matrices. In Ref. 2, p. 301, the author states that the coupler matrix is unitary. In an abstract 
mathematical sense any matrix which satisfies 

[A] (A1”T = (I] (1-A) 

is unitary when (*) implies “complex conjugate” and (T) implies “transpose.” See Ref. 6, p. 22‘). 

It can be seen that the coupler transfer matrix of equation (4) is unitary since 


C, jC^ 


K, 




jc, 


-•T 


^1 


= m 


(2-A) 


The notion of “reactive” or “lossless” is sometimes associated with the unitary matrix of a 
microwave junction, and reciprocity is evidently not a consideration. See Ref. 7, p. 273. It will be 
recalled that the objectives of the present document were restricted to “unilateral” transfer 
matrices, with the possible exception of energy conservation verification for the mismatched 
coupler. 

It is interesting to test the transmission line transfer matrix of equation (18) to detennine 
whether or not it represents a univary matrix for certain special tenninations after considering 
the general tennination. Furthermore, it appears that the “unitary” test may include “reflectionless’ 
as well as “lossless” for a given junction. 

In Ref. 8, p. 7, it is stated that the operation of taking conjugates is distributive with respect 
to addition, subtraction, multiplication, ind addition. Relying on this, the condition 




1 +ro, e 


,,p 

' ‘ 0 I 


= 1 


(3-A) 


would have to be met under a, = 0, and a similar expression for line (2) would result fora, = 0. 


A-1 


An impedance match reajlts in a zero reflection coefflcient (Pq , = 0) so that 

equation (3-A) is satisfied upon introducing a resistor. Conversely, a perfect short-circuit (Z^ = 0) 
results in a non-zero reflection coefficient (Pg| ~ -1) so that equation (3-A) is not satisfied upon 
introducing a lossless termination. 

Matrix multiplication is valid under an associative law, raising the possibility of forming 
equivalent matrices, and inquiry about the unitary character of the result. Combining the two 
square matrices of equation (-3) into an equivalent unilateral matrix is admissible, but the result 
will not be a unitary matrix. The introduction offP^j =i^0) and (P^ ^0) into equation (33) is 
an admission that energy will Ik reflected at the output end (K) of the coupler. I he discussion on 
transmission lines and unitary matrices now becomes relevant here as reflection was associated with 
a non-unitar> matrix (Z^ - 0, P^ = -1) previously. 

.Another example of a non-unitary equivalent matrix can be generated even when the reflec- 
tion coefficients at the coupler output ports are equal to zero (P,|= P ,2 “ ^)- ati attenuation 
matrix such as that of equation (30) is placed between the generator and the coupler. 
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which is not unitary due to power dissipation within the combination. 


APPENDIX B 

LOSSES IN MISMATCHED LINES 
When = Zq, equation (14) or (17) leads to 


dB = 20 lo| 


= 20 loi 




= 20 log e“® 




(l-B) 


since 


|e /^| = 

When ^ Zq, and standing waves exist on the transmission line, 


as given by equation (32). Then 


a-* a 


dB = 20 log e“^ -*■ 20 log e”' 
10 10 


(2-B) 


(3-B) 


a t 


(4-B) 


A convenient relationship may be derived relating decibels and nepers. In an abstract sense, if 

(5-B) 


and 


log N = M 

e 


log N = P 
10 


then 


e” = loP 

Taking the logarithm to the base ten of both sides, 

log N lo^ e = lo^ N . 

It follows that 


(6-B) 


(7-B) 


(8-B) 


dB = 20 log e“® “*> 20 (.434) log = 8.686 ai 8.686 a C 


(9-B) 


when ^ . 


B-1 



APPENDIX C 


CHARACTERISTIC IMPEDANCE OF LINES 


The theoiy underiying transmission lines can be found in many standard texts. See Ref. 3 and 
S. When series capacitance and shunt inductance are neglected, which is commonly done, the 
characteristic impedance is written us 


Z 


0 



(l-C) 


When the series resistance (R) is small with respect to (j^L), and the shunt conductance (G) is 
small with respect to O^^’L), the above is often simplitled to 



which implies a purely resistive characteristic impedance wh'‘n 

R = G = 0 . 


(2-G) 


(3-C) 


The propagation factor ( 7 ) is usually written as 


7 = [(R + jwL) (G +jcjC)J 




which simplifies to 


7 = a + j^ = 0+ jcu (LC) 


*/2 


(4-C) 


(5-C) 


when equation (3-C) holds. For lines with small losses, a binomial expansion (Ref. 3, p. 20) leads to 

(6-C) 


R / G /T 


In practice, a purely resistive (Z^) and a non-zero (a) are assumed even though this may appear 
paradoxical. 
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